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TOWARDS RIGOROUS ANALYSIS OF THE 
LEVITOV-MIRLIN-EVERS RECURSION 


Y. V. FYODOROV, A. KUPIAINEN, AND C. WEBB 

Abstract. This paper aims to develop a rigorous asymptotic analysis 
of an approximate renormalization group recursion for inverse participa¬ 
tion ratios P q of critical powerlaw random band matrices. The recursion 
goes back to the work by Mirlin and Evers [33 and earlier works by Lev¬ 
itov |32l [33] and is aimed to describe the ensuing multifractality of the 
eigenvectors of such matrices. We point out both similarities and dissim¬ 
ilarities between the LME recursion and those appearing in the theory of 
multiplicative cascades and branching random walks and show that the 
methods developed in those fields can be adapted to the present case. 
In particular the LME recursion is shown to exhibit a phase transition, 
which we expect is a freezing transition, where the role of temperature 
is played by the exponent q. However, the LME recursion has features 
that make its rigorous analysis considerably harder and we point out 
several open problems for further study. 


1. Introduction 

The goal of this paper is to study the asymptotic behavior of a stochastic 
branching recursion (which we refer to as the Levitov-Mirlin-Evers (LME) 
recursion) conjectured to be related to multifractal properties of the An¬ 
derson transition. More precisely, building on an earlier work of Levitov 
132 . 33] Mirlin and Evers argued in m that in the so-called critical power 
law random banded matrix model (PRBM), for small values of an auxiliary 
variance parameter, the inverse participation ratio (essentially the t v norm) 
of a typical eigenvector of the random matrix approximately satisfies the 
LME-recursion. By analyzing this recursion they further argued that the 
inverse participation ratio scales in an anomalous (multifractal) way with 
the size of the matrix. Let us also mention that a related, though less ex¬ 
plicitly developed renormalization procedure was outlined for essentially the 
same model in m, with similar conclusions. 

Our goal is to rigorously study this recursion and prove the conjectured 
scaling behavior. Our approach is based on the fact that this recursion re¬ 
sembles ones appearing in the study of weighted branching processes such 
as the branching random walk. These recursions, their applications and 
related models appear for example in other areas of probability, the study 
of algorithms, mathematical physics and finance and have been studied ex¬ 
tensively While being similar to recursions arising in 

other applications, the LME-recursion has also important differences. In 
particular, we are not able to prove in full generality the type of results 
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that are available for more traditional branching recursion and this leaves 
various open questions that might be interesting for understanding better 
the multifractal geometry of these eigenvectors. 

Thus in addition to studying the model itself, this note aims to bridge 
the gap between the Theoretical Physics and Probability communities in 
both directions. Namely, on one hand we are presenting to readers from the 
physics community with some of the results and tools used in the Probability 
communities to study some multifractal geometric objects. On the other 
hand we present to readers from the Mathematics community with a new 
type of stochastic recursion which has important applications in random 
matrix theory and its relations to physics of disordered systems and is an 
interesting and rich object in itself. In particular, there are many open 
questions we are not yet able to answer. 

In Section 2 we briefly recall the critical power law random band matrix 
model and the heuristic derivation of the LME recursion for the inverse 
participation ratios. In Section 3 we recall branching recursions arising in 
the study of branching random walks and the type of results and tools one 
has there. Section 4 studies the scaling properties of the moments of the 
solution of the LME recursion. We deduce the existence of a phase transition 
in the model, namely for small values of the index of the participation ratio, 
the participation ratio lives on the scale of its mean, while for large values 
it does not. Using the knowledge on the scaling behavior of the moments, 
we are able to prove (using the method of moments) that for small values of 
the participation ratio index, the correctly normalized inverse participation 
ratio converges in law to a non-trivial random variable whose law can be 
characterized in terms of a non-linear integro-differential equation. Finally, 
in Section 8 we discuss some open questions suggested by the connection 
with branching random walks. 
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2. The model 


To motivate our stochastic recursion, we recall the model and argument 
presented in m, based on the earlier work in [32]. We consider a power-law 
random band matrix model [38] - that is random N X N real symmetric ma¬ 
trices H with independent entries whose variance decays with the distance 
from the diagonal in a power law fashion: 


( 1 ) 
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with b being a coupling constant controlling for a fixed a the relative strength 
of fluctuations of the entries close to the main diagonal and away from it. 
In particular, for 5 > 1 there is a well-defined band of the width b around 
the diagonal where the variance remains constant, and then decays in a 
powerlaw way outside that band. In the opposite case b —> 0 the condition 
| i — j | < b implies i = j so that the matrix has a dominating diagonal with a 
small admixture of off-diagonal entries with powerlaw-decaying amplitude. 

Such a model for a = 1 is considered in Theoretical Physics literature as 
a paradigmatic example of critical random matrix models with multifractal 
eigenvectors, see e.g. [31] for a recent discussion. As such it played in the 
last decade an important role in attempts to understand generic proper¬ 
ties of wave functions for systems at the point of the Anderson localization 
transition where such multifractality is conjectured to be one of the most 
characteristic features m- Let us note that a few other critical random 
matrix models with multifractal eigenvectors attracted much of attention 
recently mmm- 

It was first suggested in [38] that in the limit of large control parame¬ 
ter 1) > 1 and Gaussian randomness the random matrix model 0 under¬ 
goes an Anderson-type transition when changing the exponent parameter a. 
Namely, via a mapping to a non-linear a-model it was argued that whereas 
for any a > 1 the eigenvectors of the underlying matrices are localized, they 
are delocalized for a < 1, and precisely at the critical value a = 1 eigen¬ 
vectors show nontrivial multifractal behavior. To quantify this statement 
consider a typical eigenvector -i/’ = (ifii, ..., i/t/v) and given a parameter q > 0 
define the inverse participation ratio in terms of the l 2q norm as 

N 

( 2 ) P(q,iP) = H\\l g q = ^2\H 2g - 

i =1 

where we normalize the norm to 1 i.e. P(l,ip) = 1. One expects P(q) to 
scale with volume as 

(3) P(q,^) oc AT _d O)(<?- 1 ) 

with d(q ) = 0 for localized states, d(q) = 1 for extended states and 0 < 
d(q) < 1 for multifractal states. 

We consider in this paper only the critical case a = 1 in 0 leading to 
multifractal eigenfunctions. Let us also for simplicity take periodic boundary 
conditions i.e. we let the indices j E Zjv, the integers modulo N, and let 
| i — j | in 0 stand for the distance on Zjy and similarly the addition i + j . 
The LME renormalization group deals with the same model ([I]) in the case of 
a small control parameter b -A 0 which is exactly opposite to that considered 
in [38]. Nevertheless, by heuristically deriving and analyzing an approximate 
renormalization group flow Mirlin and Evers concluded that for the critical 
value a = la nontrivial multifractality of eigenvectors survives for any small 
b> 0, although all d{q) in this limit remain, non-surprisingly, parametrically 
close to the localized limit d(q) = 0. 

Essentially, the LME procedure aims to construct the orthogonal matrix 
O diagonalizing H, O r HO = diag (Ei) inductively by finding matrices O n 
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S.t. 

(O^HO n )ij = 0, 1 < \i - j\ < n. 

Obviously 0$ = 1 so let n = 1. We have that Ej = Ha are i.i.d with 
variance 1. Hence typical eigenvalue differences are E'j — Ej = 0(1). More¬ 
over, the variance of the perturbation hij := Hijlu_ji = i is b 2 so that typi¬ 
cally hij = 0(b). If these conditions were satisfied the matrix 0\ would be 
0\ = l-|-m where m = 0(b). Indeed, to first order in the perturbation h 

( 4 ) = E o h ^ E o + • • • • 

3 1 

However, when Ej — Ej = 0(b) perturbation theory breaks down, and such 
an event is called a resonance. Pick a > 0 and define the resonance set at 
scale 1 by 

R 1 = {i:3j,\i-j\ = l,\Ej-Ej\<b a }. 

With high probability, as b —> 0 the set i?i consists of isolated nearest 
neighbor pairs separated by distance 0(b~ a ). Under such an event the 
matrix h is diagonalized by the matrix 

(5) 0 Rl := 


where 0\ diagonalizes the resonant matrix 

( H ii+ 1 \ 

V E l i ) ' 


The LME RG will now be defined by making the approximation 


(7) O i = 0 Rl © l R c 

and approximating OjHO \ by diag (Ej) + H 1 where Ej = Ej when i G Rq 
and Ej, Ej +1 are the eigenvalues of ([b]) if i, i + 1 E R\ and Hj- = 

Note that the true O i has of course a more complicated structure. In 
the non resonant region its off-diagonal matrix elements O \ E will not be 
exactly zero, but small and exponentially decaying in the separation | i — 
j |. Moreover the resonant and non resonant regions are not completely 
decoupled. However for small b Q should capture the essence of 0\ . Also, 
OjHO i will include matrix elements with \i—j\ = 1 and further renormalize 
the diagonal elements Ej. These corrections are subleading in b and are 
dropped in the LME approximation. 

Inductively, we approximate Oj(HO n by diag(E’”) + H n where Hj- = 
Hijlu_j\> n . We define the resonance set at scale n+1 by 

R n+ i = {i : 3j, \i-j\=n + 1, | Ej - Ej\ < (^)“ }. 


and set 


O n +i — O n o n 


where o n is defined as 0\ above by diagonalizing on the resonant set R n +i 
the matrix 


M n . ■ = 
13 


Ej hij \ 
hij Ej ) 


(8) 
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where hij = -^T)l|i-j|=n+i • hij have variance 6 2 /(n + l) 2 . This time the 
resonant set has density 0((^pp) a ). 

If we set H n = 0 the eigenvectors of H are given by ip™ = O n ip® where ip® 
are the canonical basis vectors of M^. In case of a resonant pair i , j write 
the restriction of o n to V vj := span{^,^°} as 


(9) 

so that 


°n\Vij 


cos 9 n sin 6 n 
— sin 9 n cos 9 n 


ipi +1 = cos 9 n tpi + sin 9 n tp'j 
ip™ +1 = cos 9 n ip 1 j - sin 9 n ipi. 

We point out that the random variable 9 n here depends on the resonant 
pair (i,j), but for all pairs, the distribution is the same. 

Hence we get for the the inverse participation ratios 

P(q\ € +1 ) = II cos 9 n ipf + sm9 n ipy\\ 2 2 q q 

Note that ip™ has support in the original basis in a ball of radius k around 
i where k is the largest scale k < n s.t. 'iff was resonant. Hence with high 
probability ipf and ip™ have disjoint supports and on that event 

(10) P(q;W +1 ) = (cos 2 9 n yP(q-W) + (sin 2 9 n fP{q^). 


Remark 1. The terms (sin 2 9 n ) q and (cos 2 9 n ) q will appear repeatedly and 
for notational brevity, we will write them as sin 2 * 7 9 n and cos 2q 9 n (the possible 
ambiguity in the latter notation being when sin 9 n or cos 9 n are negative). 


The eigenvalues Ef and E™ are functions of the original EHi m with 
disjoint support in the indices and hence independent in our approximation. 
Obviously the law of P(q;ip™) is independent of i. We are therefore led to 
consider the iteration for random variables P n (q)- 


(ID 


Pn+i(q) = sin2<? QnP^Xq) + COS 2q 9 n P^\q), 


where Pn^ are mutually independent copies of P n {q ), independent of 9 n and 
Pi(q) = 1. Note that (11) differs from (10) in that the latter was argued 
to hold (approximately) provided the resonant set consists of disjoint pairs 
whereas in the former we extend that relation everywhere. In fact shortly 
we will argue that this change does not affect the multifractal behavior. 

To understand the recursion ( [TT| ), let us derive the law of 9 n . First note 
that 9 n , 9 n + bn and — 9 n lead to the same eigenspaces. Hence we may 
restrict 9 n to the interval 9 n G [—f, f]- Introducing the variable 

Tn = 5E n /hij 

where 5E n = b{Ef — E '") the eigenvectors v = (a, b) T of ([8]) satisfy 
(r n ± sgn(h^)v / l + r 2 )a + 6 = 0. 


T n = 


Writing 

( 12 ) 


cot 29 
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some algebra leads to | = — cot 0 or | = tan 0. Hence we may identify 6 
with 0 n . 

Let the law of t := 5E n have density p{t ) and let h tJ have the law of ev 
where e = ^ and v has density j(v). We suppose p and 7 have all positive 
moments and that H/O^oo < 00 . Then for a bounded continuous function / 

(13) E(f(9 n )) = f dtp(t ) f dir/(v)f (-^cot _1 (^)) 

Jr Jr ' 7 

Changing variables by v = — | tan 29 we get 


(14) 

where 

(15) 


E(/(0n)) 


2 /■? d0 
e J_k cos 2 20 

4 





ri0 

sin 2 29 


Xe(9)f(9) 


tan2l f )f(9) 


Xe{9) = 2 f dtp(-^gt)\t\j(t). 
J M 


Note that Xe regularizes the integrand in © around the singularity at 
0 = 0. Indeed, by change of variables 


(16) X.(D) < 2(S^) 2 ||7lU f dt p (t)\t\ < 

J R 


Moreover, we have 


Lemma 2. Let f be bounded on [—f, f] with f(9) = 0(\9\ a ) as 9 —>• 0, with 
a > 1. T/ien 

E(/(*n))=^ / 4 T^/(0) + O((!r) 

n ,/_i sin 2t/ n 

4 

with (3 = 2 p( 0 ) J R di|i| 7 (t). 

Proof. We have by © 

6 / < Ce- 1 ^ \9\ a < Ce a . 


Moreover, by our assumption on /, 

f d9 


sin 2 20 


/(0)l|tan20|<e < Ce° 


From (15) we have 
(17) 

so that 


\Xe{8) ~P\< C^^Wp'Woo < C||| 


d0 


£ J ^29 lXe{9) ~ P\f^ l \^e\>e < Ce a . 

Putting things together gives the claim. 


□ 
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Remark 3. The parameter j3 depends on p(0) i.e. the density of the eigen¬ 
value differences 5E n at 5E n = 0 and on E|u|. Note that we have (in our 
approximation) 

\5E n+1 \ = y/6E* + hl 

where h n = hij with \i — j\ = n. Therefore f3 n tends to a limit f3 as n —> oo. 
The multifractal exponents turn out to be proportional to (5. Since it plays 
no further role in the analysis we will set it to one in what follows. This 
will change some of our results from those in m by a factor of 8/ir. 

Motivated by these considerations we can now define the recursion which 
we will study in this paper. 

Definition 4. We say that the law of P n (q) satisfies the LME recursion if 

(18) Pn+i(q) = sin 29 e n PP(q) + cos 29 9 n P( 2 \q), 

where P^ and Pn^ are independent copies of P n (q), Pi(q) = 1 and 9 n is 
independent of Pn\q) and distributed on [— with density 

(19) r n (9) = (0)(sin 2#)~ 2 

and x is any non-negative function s.t. f r n = 1 and ( flfi] ) and © hold. 
We also recall that by sin 29 9 n and cos 2q 9 n we mean (sin 2 9 n ) q and (cos 2 9 n ) q . 

Remark 5. Note that as n —> oo, the law of 9 n degenerates to 9 = 0 a.s. 

Remark 6. If we write f n for the density of the law of P n , then the LME 
recursion in terms of f n is 

/•oo /♦oo 

d9r n {9) / / dpidp 2 fn{pi)fn{jp 2 )d{p - pis\VL 2q 9 -p 2 cos 2q 

Jo Jo 

which is essentially a discrete version of equation (43) in [37]. Similarly we 
can express the recursion in terms of the Laplace transform of the law of P n . 
Let us write ip n {t) = E(e~ tPn ). Plugging the recursion into this gives 

/ n 

r n {9)if n {t sin 29 9)if> n (t cos 29 9)d9. 

-f 

Remark 7. Recall that we defined the resonant pairs by \Ef — E~j\ < ( b/n) a . 
Since hij = (D(l/n) this implies that for non-resonant pairs \9 n \ = 0{n a ~ l ) 
and therefore provided 2q > 1 which we will assume throughout the paper 
we may arrange 2q(l — a) > 1 and 2(1 — a) > 1. This means that the 
contribution of the non-resonant pairs in the iteration © is 0{n~ l ~ a 'P n ) 
with a' > 0 which will not contribute to the multifractal exponent. 

Remark 8. Our approximations in deriving the LME RG were dropping 
corrections to the renormalized eigenvalues Ef +1 and the matrix o n in the 
non resonant region and ignoring multiple resonances. The perturbation 
theory in the non resonant region produces off-diagonal terms to o n that 
are exponentially decaying (at the scale n) and down by inverse powers in n. 
Similarly, Ef +1 receive corrections that are non-local in the El and h. Thus, 
even if no resonances were present the eigenvalues Ef are not independent 
but weakly correlated. In the presence of resonances perturbation theory in 


fn+\{p) = f 
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the non resonant region does not decouple from the resonant one and there 
are further corrections that should be small due to the small probability of 
resonances. However to make a rigorous proof is very challenging. It would 
be very interesting to try to adopt the iterative scheme ofm to this problem. 

3. Some background in recursions of branching random walks 

Let us consider an extremely simple example of a branching recursion 
related to a branching random walk. Let Vj and V 2 be i.i.d. standard normal 
random variables (compared to these play the role of log sin 2 9 n and 
log cos 2 6 n ), = 1 (Z corresponding to P), (3 E M (corresponding to q) 

and define 


(21) Z n+1 (P) = e«2W(/i)+e«^(/3), 

where ' are independent copies of Z n and they are independent from the 
Vi as well. Compared to , we point out immediately some differences. In 
the terms corresponding to (e^ Vl ,e^ V2 ) are not independent, moreover 
their law depends on n and as n —> 00 and their joint law converges to 
a degenrate distribution (one of them is almost surely one and the other 
almost surely zero). 


The asymptotic behavior of Z n and related objects has been studied under 
various generalizations of this recursion, see e.g. 13123 EDI Q21 M ED El 113 
E3 03 E3IIS103]. Gaussianity is not important, one might have a random 
amount of terms in the recursion and the terms might even be correlated. 
Though much of the behavior in such recursions is universal (21) is perhaps 
the simplest recursion to study these universal features in. Moreover, even 
situations where the law of the terms corresponding to V) depends on n has 
been studied in some relatively simple cases, but it seems that the situation 
where the limiting law of these has such a degenerate structure has not been 
studied rigorously. 


Before briefly discussing the type of methods used to study recursions of 
this type, let us state the theorem describing the asymptotic behavior of 
Zn(P). 


Theorem 9 ( [13 @ US EJ ESI 03 E] ). For /3 < (3 C = v / 2Toj ? 2, M n (/3) = 

E = 2 -n e -n ^ Z n (f3) converges in law to an almost surely positive 
random variable M((3), which satisfies the distributional equation 


(22) M(/3) = + i e ^-irM (2 )(/3), 

where are independent copies of M(fj) and are independent ofVj. 

For (3 > Pc, M n (P) converges in law to zero, but yHiM n (f3 c ) converges in 
law to an almost surely positive random variable M'(f3 c ) which satisfies 


(23) M'(/3) = \e ^ Vl ( M ' (/3))( 2 ), 
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with similar independence structure. 

Finally, for (3 > /3 C , M n (/3) converges in law to a non-trivial 

random variable M(f3) satisfying 

(24) M(0) = I e ^i-^+%M( 1 )(/3) + ^e 0V2 ~^ c+ ^M {2) (/3). 

Moreover, the law of M ((3) can be written as 

(25) M(/3) = L & (M'(/3 C )), 

|8 

where the process t ha L a (t), t > 0 is a stable Levy subordinator of index a, 
independent of M'((3 C ). 

Remark 10. As we shall soon see, Z n $ can be viewed as the partition 
function of a type of Random Energy Model with specially (logarithmically) 
correlated energies and the fact that the law of M(/3) is essentially charac¬ 
terized by M’(f3 c ) is one aspect of the freezing transition occurring in the 
study of log- correlated disordered systems, see mmm for discussions of 
various physical aspects of freezing and its relation to multifractality, and 
further references. 

Remark 11. For this theorem, the relevant property of the stable subor¬ 
dinator is that it satisfies K(e~ sLa ^) = e~ ctsa for some positive constant 

C. 

In addition to proving convergence, the following properties of the limiting 
random variables are known: 

Theorem 12 ([27, 13], [T8J |39j). For [3 < (3 C , E(M(/3) P ) < oo if and only 

o2 

tf P < pi and IE(M'(/3 c ) p ) < oo if and only if p < 1. In fact, one has the 
sharper result that there exists a positive constant c(/3) so that for (3 < f3 c 


(26) 

£ 

lim xP 2 P(M(/3) > x 

) = C(/ 3) 

£—>•00 

and 

(27) 

lim xP(M'(/3 c ) > x) 

x —>oo 

/ ""1J 

TL 

II 


— _& 

Remark 13. For (3 > (3 C , the tail asymptotics P (M(/3) > x) ~ x P fol¬ 
low from the representation M((3 ) = Lp^{M\(3 c )) and the corresponding 

p 

asymptotics of L ^(1) and M'((3 C ). 


Let us now turn to a few words on the tools typically used for studying the 
recursion. Perhaps the most efficient (and most easily generalized) method 
of studying Z n ((3 ) is through a branching random walk representation. Con¬ 
sider a binary tree and on each edge, place i.i.d. copies of standard normal 
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Figure 1. A graphical representation of the binary tree and 
the related random variables used to construct a branching 
random walk. Here Aooi = Vo + Voo + Hjoi and Adoi = 
V\ + Vio + kioi- 


variables. Each edge at hight n in the tree can be labelled by a binary 
sequence a = (cti,..., a n ) G {0, l} n , where 0 refers to taking the left path 
and 1 to taking the right one at a branching point. Let us denote by V a the 
variable associated to the edge a and if a G {0, l} n and k < n, write a\k for 
the truncation of a to length k: a\k = (ui,..., ct/ c ). 

Let us truncate the tree at height to. Each branch in the truncated tree can 
again be identified with a binary sequence a = (<7i,..., a n ) G {0, l} n . Call X a 
the sum of the Vj along the branch: X a = ^o-l k- The collection (X a ) 

can be viewed as a branching random walk. See Figure 1 for an illustration 
of the construction. Now let us define 


(28) 


MP)= E 

o-ejo,!}" 




Remark 14. Note that if we think of —X a as the (random) energy of a spin- 
configuration a, Z n {j3) becomes the partition function of a kind of Random 
Energy Model with correlations. In fact, one can check that the covariance of 
X a and X a t is the logarithm of the distance of a and a' in the dyadic metric 
of the tree. Partition functions of Random Energy models with logarithmic 
correlations in Euclidean metric have attracted recently considerable atten¬ 
tion due to several interesting applications in Statistical Mechanics, Random 
Matrix Theory and Number Theory see e.g. ED E2 (231251E] and references 
therein. They share many characteristic properties with their diadic tree 


counterparts, such e.g. as the tail behaviour of moments (26). Analysis of 


such models frequently proceeds by identifying an approximate binary tree 
structure, see e.g. [3]- 

For a G {0, l} fc , we also define a similar object corresponding to the tree 
of height n rooted at a: 
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(29) = Y 

cr'e{ 0,l} n+k - 

<j'\k=(T 

Note that Z^\[3) = Z n ((3) and for a, a’ € {0,l} fc and er ^ a', Z^\f3) 
is independent from Zi a \/3) and X a is independent from Z^ 7 \/3) for any 
a, a' 6 {0, l} fc . Also one has (simply plugging in the definitions) 

(30) Z n+k (P)= Y e? Xa ZP{P). 

o-e{o,i} fe 

In particular, specializing to k = 1, 


(31) 


J n +1 


(/3) = e^Z^(/3) + e^Z^(i3), 


so we have realized the variables appearing in our recursion on the same 
probability space. Moreover, if we denote by T n the u-algebra generated by 
the V a with a G {0,l} fc for any k < n, then it follows directly from (30) 
(by setting n = 1) that M n ((3) = a positive martingale and one 

sees immediately that it converges to some limit. Then analyzing moments 
of M n (/3) through the recursion, one finds that the limit is non-trivial for 
(3 < (3 C and zero for /3 > j3 c - Our analysis of the moments of P n (q) will be 
very similar to the way one can analyze the moments here so we shall not 
discuss this further now. 

We do however point out that while in the model we have described, the 
limiting object M(/3) has moments only up to order p < /3^//3 2 , Barral and 
Mandelbrot [Bj considered a slightly different kind of model where e^ Vi were 
replaced by variables that could take negative values as well, and in this case 
one has a limiting random variable which has moments of all positive integer 
order. Moreover, they used the values of the limiting moments to identify 
the law of the limiting random variable. This is similar to our argument for 
small q. 


Finding the correct normalization for (3 > f3 c is far more complicated. 
At f3 c what turns out to be critical is proving that the so-called derivative 
martingale (see jlO] ) 


(32) -M' n (p c ) = 2~ n Y (nf3 c -X a y* x °^ n 

a£{ 0 , 1 }" 

converges to an almost surely positive random variable. In fact, the limit 
equals (in distribution) M'(/3 C ) up to a constant multiple. 

For /3 > /? c , what turns out to be the important question is the behavior of 
maXo-gjoqjn X a . What is typical of the freezing transition is that for large /3, 
only the variables X a which are close to the maximum one (on that level in 
the tree) contribute to Z n (j3). The result for the maximum is the following: 
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Theorem 15 (fl* S3). There exists a positive constant c > 0 such that 


(33) 


lim P 

n—foe 


max X a — v/2 log 2 n H- . 

ae{o,i}™ V 2y / 21og2 


logn < x 


= E(exp (-ce-^ x M'{/3 c ))). 


Remark 16. One feature of this limiting distribution is that it has a tail 
of the formm xe~^ cX (as x -A oo) so this is a problem of extreme value 
statistics that is not in the Gumbel universality class. As conjectured in the 
physical literature [HI [22] such a tail together with the coefficient 3/2 in (33) 


are characteristic features of the generic sequences, processes and fields with 
logarithmic correlations, and are intimately related to the powerlaw tail of 
the partition function moments (26), see [23]. Rigorous results corroborating 
this conjecture can be found e.g. in [16] and references therein. 


The proof of Theorem 15 is quite technical, but a fundamental idea is 
proving that branches X a that live on the scale of the maximum end up 
there along paths which perform an excursion around a linear path of speed 

\/2bp. 


4. Scaling of moments in the LME recursion 


Compared to the branching random walk, there is an immediate difficulty 


in analyzing the convergence of the solution to the LME recursion (18). 


Namely even if one realizes the random variables on the same probability 
space in tree form, one does not have a martingale structure. We have 
not noticed any monotonicity or really any contractivity either that would 
enable proving convergence in general. Thus instead of focusing on proving 
convergence, we will begin by studying the moments of P n (q) satisfying the 
LME recursion (18). Using arguments similar to ones used for the branching 


random walk, we will see that there is a phase transition: there is an explicit 
q c (approximately 2.4056) such that for | < q < q c , for suitable p > 1, 

Pn(q) 


E (P n (q) p ) ~ E (Pn(q)) p While for q > q c , 


E(Pn(<?)) 


0. In the following 


sections we will also use these results to show that for \ < q < 1, 


2 " * E (P n (q)) 

converges to a non-trivial random variable. For brevity, let us introduce a 
name for this random variable. 


Definition 17. Let P n (q ) be given by the LME-recursion (18) and let q > 
We then write 

Pn(q) 


(34) 


n n(q) = 


HPn(q))' 


Note that if one is able to prove that sup n E(II n (g) p ) < oo for some p > 1, 
then one has immediately that there exist constants 0 < Ci, C 2 < 1 such 
that for all n, 


(35) 


P(n n (q) > 2 ) < Cl and P (n n (g) > 5 ) > C 2 . 















ANALYSIS OF THE LME RECURSION 


13 


The first inequality is simply Markov’s inequality, since E(n n (g)) = 1. 
The second inequality is a Paley-Zygmund type inequality: we have by 
Holder’s inequality 


1 = E ^n„(g)l jn n (g) > + E ^n n (g)l jn n (<?) < 

<(E(n n (gf))^p(n n (g)>l)y P + 1 


1 

2 


yielding 


(36) 


n„(g) > - ) > 


1 


2(E(IL n (q)P))p 


V 

P-1 


So in particular, the scaling properties of the moments imply that we can’t 
have n n (g) ~ oo with high probability or n n (g) ~ 0 with high probability. 
Thus for such q , E (P n (q)) is a scale (though we expect it’s the only one) at 
which P n (q) lives. 


For our proof, it will be important to understand the asymptotics of 
E (P n (q)). To do this, we consider T n (q) := E ^p +1 ^' 1 = E(sin 2 " # n +cos 2 " # n ). 
The central result for this quantity is 

Lemma 18. For q> we have for a = min(2, 2g) 


(37) 1 - T n {q) = E(1 - sin 2 " 6 n - cos 2 " B n ) 

where 


b T(q ) + O 
n 



(38) 


T(q) 


8 T(g — 1) 


Remark 19. Note that T(q) is finite for q > 


Remark 20. We choose to write the errors in terms of - instead of - as 

J n J n 

one might potentially be interested in making b depend on n. 


Proof. Note that as 0 —> 0, 1 — sin 2 " 0 — cos 2 " 6 = 0(|#| mm ( 2 "’ 2 )), so by 
Lemma [2j 


7 r 

(39) 1-T n (q) = - [ 4 - 7 -^—(1 — sin 2 " 9 — cos 2 " 9) + 0((fm -1 )“). 

n J_iL sir 20 
4 

Recalling our convention for sin 2 " 6 and cos 2 " 0, we see that the integrand 
is symmetric under 6 —> —9. We also note that for 9 6 [0,7 t/4 ], (1 — cos 2 "0 — 
sin 2 " 9)/ sin 2 29 is symmetric under 9 —> —9 + ir/2 so we find that 


(40) 


T <«) = 4 


-P 

4 to 


cW 


sin 2 9 cos 2 9 


(1 - sin 2q 9- 


cos 


2 q 


9). 
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To see that the integral can indeed be written as a ratio of T-functions, 
let us introduce an auxiliary complex parameter to the integral: 


(41) T(q, z) := - f (1 — sin 29 0 — cos 2<? 0)(sin 0 cos 0)“ 2 d0. 

4 Jo 

We see that T(q, 0) = T(q) and T(q,z ) is analytic in z for Re(z) > 
max(l — 2g, —1). For Re(z) large enough, we see that each term in the 
integral can be expressed as a Beta-function, which we can write as a ratio 
of T-functions: for Re(z) > 1, 


(42) 


_ i r(^) 2 _ i r(g + ^)r(^) 

8r(;s-l) 4 T(q + z-l) 


As the only singularities of the T-function in the complex plane are simple 
poles at non-positive integers, we see that for Re(z) > max(l — 2q, — 1), the 
only possible singularity of the meromorphic function 


i r(f) 2 i r( g + ^)r( 

8r(z-i) 4 r (<7 + z- 

is at z = 1, but we see that at z = 1, the poles of the two different terms 
cancel and 




(44) T(q,z) = 

even for Re(z) > max(l — 


i r(f) 2 i r( g + ^)r(^i) 

8T(z-l) 4 T{q + z-l) 

2q, —1) . In particular, 


(45) 


T(q) = T(q, 0) 


i r( g -i)r(-i) y^ T(q-l) 

4 T(q-l) 8 T(q- 1)' 


□ 


We shall also need similar estimates for joint moments of sin 9 n and cos 9 n . 
Lemma 21. For q> \ and j, k > 1 we have for a = min(2 qj, 2 qk) > 1 

(46) E(sin 29 - J 0 n cos 2qk 9 n ) = — ^ J sin 2q ^ 2 9 cos 2gk ~ 2 9d9 + O 

Proof. Noting that as 9 —> 0, sin 29 - 7 0 cos 2<?fc 0 = 0(\9\ a ), we again refer to 
Lemma [2] We also make use of the symmetry under 0 —0 (hidden in our 

notational convention). □ 


Before moving on to studying the moments of P n , we fix some notation. 

Lemma 22. There is only a single solution (which we call q c from now on) 
to the equation qT'(q) = T(q) in the domain (^,00). 
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Proof. As this is equivalent to showing that the function H (q) = qT'(q) — 
T(q) has only a single zero, it is enough to check that H'(q ) = qT"(q) does 
not change its sign and that H(q) does. The first claim is clear for example 
from the representation 

(47) T(q) = - / . (1 - sin 2 » 6 - cos 2 » 9). 

4 J o sin 0 cos^ 9 

For the second one we note that as T(q) < 0 for q < 1 and T'(q) > 0 for 
all q. we see that qT'(q) —T(q) > 0 for small q. For large q, one can use the 
asymptotics of the polyganuna function of order 0 to check that H(q) < 0 
for large q (in fact, one can check that H(q) ~ — \\[q for large q). Moreover, 
numerically one finds that the root q c is at 2.4056. □ 

We also need yet another elementary remark to describe which moments 
of n n (g) are bounded 

Lemma 23. For \ < q < q c , the set {p > 1 : T{pq) —pT(q) > 0} is an open 
interval. 

Proof. For g(p) := T(pq) - pT(q), g'(p ) = qT'(pq) - T{q) and g"(jp) = 
q 2 T"{pq). If the set was not an interval, there would be several points at 
which g'(p) = 0. But since g" is strictly negative, g' is strictly decreasing, 
we see that this is impossible. The interval is non-empty since for q < q c , 
</(!) > 0 . □ 

Remark 24. For q < q c , we will write {p > 1 : T(pq ) — pT(q) > 0} as 
(1 ,p*(q))- Moreover, as d q (T(pq)—pT(q)) = p(T’(pq) — T'(q)) < 0 for p > 1 
(as T"(q) < 0), we see that p*(q ) is decreasing in q. 

We are now in a position to characterize which moments are finite: 

Proposition 25. For \ < q < q c , sup n E(II n (q) p ) < oo for 1 < p < p*{q). 
For p > p*{q), E(n n (g) p ) —> oo as n -* oo. 

Proof. Let us begin by considering p < p*(q). By Jensen’s inequality, for 

p' < p, E(n„(( 7 ) p ) = E(n n ((/) p p) < (E(II n (g) p )) p . Thus it is enough to 

consider p close enough to p*{q). Let us further assume that p < 2 for now 

£ 

and use the subadditivity of x e->- x 2 to see that 


E (P n+1 (q) p ) = E ((sin^pW^) + cos 2 * 0 n pW (q)fi) 

< E(sin 2 ® 0 n P^ ( q) p + 2 sin® 0 n cos® 0 n P.^ (q) § P^ (q) 2 ) 

+ E(cos 2 ® 0 n Pj l 2 \q) p ). 

Normalizing by K(P n+ i(q)) p and noting that by Jensen, E(P n (q) 2 ) < 

P 

E {P n (q))z, we see that 


(48) 

E(U n+1 (q) p ) < 


E(sin 2 ® 0 n + cos 2 ® 0 n ) E(2 sin® 9 n cos® 0 n ) 

(E(sin 2? 6 n + cos 2q 9 n )) p " ^ (E(sin 2<? 6 n + cos 2q 9 n )) p 
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Let us now assume that E(II n+ i((/) p ) > E(n n (g) p ). We then obtain that 


(49) 

nf E(sin 2 ®^ n + cos 2 ®^,) \ n E(2 sin® 2 6 n cos® 2 0 n ) 

b \ (E(sin 2l? + cos 2l? # n )) p / ,l+1 ^ — b (E(sin 2<? 9 n + cos 2 * 2 9 n ))P 


As n —> oo, the term multiplying E(n n+ i(g) p ) converges to T{jpq ) — pT(q) 
which is positive by our assumption. If q > 1, the right side converges (by 
Lemma 21) to 


(50) f sin®’ 2 0cos qp 2 Odd. 

Jo 

and this integral converges for any p > 1. If q < 1, we note T(pq) — 
pT{q) > 0 for any p > 1, in particular choosing p = 2, we can again apply 
Lemma Kill and we get a convergent integral as q > We conclude that 
if E(n„ 4 li(<?) p ) > E(II n (g) p ), we get an absolute bound for E(IL rt+ i(g) p ). 
On the other hand, if E(n ra+ i(g) p ) < E(II n (g) p ), we keep going back in the 
sequence until we have an index of increase or we reach n = 1 - in any 
event, we get an upper bound. We conclude that E(II n (g) p ) is bounded for 
pe [1,2]. 

One can then proceed inductively to prove this if p £ (A: — 1, k], for k > 2: 
one writes p = k^ and with a similar subadditivity argument one has 

(51) E(n n+1 ( g ) p ) < Q^(sm 2qp ie n cos 2 ^e n ) 

xE(n n ( g ))f)E(n„,(g)^). 

Assuming E(II ri+ i(g) p ) > E(II n (g) p ) gives 


(52) 

E(n„ +1 (^) p ) (i 


T n (qp) \ 

Tn(q) p J 


s ^ E (*) E ( sin2 "‘"») 

xE(n n ( 9 ))¥)E(n n ( 9 )^). 


We now assume that we know that sup n E(n r j(g) p/ ) < oo for 1 < p' < k— 1. 
As p < k, we can again use Jensen to estimate E(II n (g)v) < E(II n ((jf) J ")£ 
which is bounded in n for 1 < j < k (by our induction hypothesis). We can 
then proceed as for p < 2. 

Consider now p > p*(q). Then by superadditivity of x H > x p for p > 1, 
we have 


E(n n+ 1 (g) p ) > 


E(sin 2 ®0 n + cos 29p 0 n ) 
(E(sin 2,? 6 n + cos 29 9 n ))P 


E(n n (g) p ). 


(53) 
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The prefactor is 1 + ^(pT(q) — T(pq))+0((^) l+e ) for some e > 0 so iterat¬ 
ing this, we find that E(II n (g) p ) > (7 e bl °g n (p T (i)- T (P9)). By our assumptions, 
the exponent is positive so we see that this blows up. 

□ 

Remark 26 (Multifractal scaling). We point out that this result implies the 
multifractal scaling we discussed earlier (see ©/• for q < q c , E(II n (q) p ) is 
bounded for some p > 1 and we argued that this implies that asymptotically, 
n n (q) is non-trivial. In fact, we expect that it converges, so we’ll make the 
assumption that as n —>• oo ; II n (q) = 0(1 ) (we’ll see that for q < 1, we have 
convergence so at least in this case, our argument is solid). We thus have 


P n (q) = U n (q)E(P n (q)) 

n / 1 \ 

(54) =0 (1 )II 1-1%) 

k=l ' ' 

= 0(l)n~ bT{q) . 

In the notation of ([3]), this implies 


(55) 


d(q) 


bT(q ) = x/tFT (q - \) 
q- 1 8 T (q) 


If we fix some q > 1/2, by taking b small enough, we have d(q) £ (0,1) 
implying the claimed multifractal scaling for P n (q)- For larger values of q, 
we expect a different normalizing factor, but still one that produces multi- 
fractality. 


To see that q c described here is truly the critical point in that for q > q c , 
normalizing by the expectation is not the correct way to proceed, we point 
out the following result: 

Proposition 27. For q > q c , there exists an h = h(q) < 1 such that 


(56) E(n n (q) h ) —1 0 

as n —>• oo. 

Proof This follows essentially from the subadditivity of x i— > x h . We get 
immediately that 


(57) 


E(U n+1 (q) h ) <E(U n (q) h ) 


E(sm 2qh 9 n + cos 2 ^ On) 
(E(sin 2l? 9 n + cos 2 ® 0 n )) h 


The ratio here is 1 — ^(T(cjh) — hT(q )) + 0((^) 1+t ) (for some e > 0). As 
qT'(q) — T(q) < 0, One can check that for h 6 (0,1) close enough to one, 
T(qh) — hT(q ) > 0. Thus 


(58) 


as n —)> oo. 


E(n n (g) /l ) < Ce- blo ^ T ^~ hT ^ —t 0 


□ 
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5. Convergence of positive integer moments 

In this section we will prove that for each positive integer k and for 
q € {\,qk), where qk is the unique point satisfying T(kqk) — kT(qk) = 0, 
E(n n (g)- J ) converges for j < k. 

Proposition 28. For any fixed positive value of b, \ < q < qk and j < k, 
E(n n (g)- J ) converges to a positive number Mj = Mj(q) and for l < k one 
has the recursion 


(59) Af = — - 1 -— V [ 4 sin 2 «'- 2 0cos 2q ^~ 2 6d9 

T{lq) ~ IT(q) “ \j) J lj 2 J 0 

with Mi(q) = 1. 

Proof. Let us begin with the k = 2 case. Normalizing the recursion ©by 
expectations, we find that the second moment satisfies the recursion 


E(n n+1 ( g ) 2 ) = E(n n (g0 2 ) + - i) E(n n (g) 2 ) 

+ rr T N2 E ( sin2g COS 2q e n ). 

T n {qr 


By Lemma 18 and Lemma [2 1 [ we can write (for some a > 1) 


(60) 


n 


-(E(n n+1 (g) 2 - E(n n (g) 2 )) = / sin 2< ?- 2 0cos 29 - 2 ^ 

b Jo 

- (T(2q) - 2T(q))E(U n (qf) + O ( 


a— 1' 


From this we can see what the limit of E(II n ((/) 2 ) must be if it exists. If 
E(II n (g) 2 ) were to converge, the left side can only converge to zero - if it had 
a non-zero limit, E(II n (g) 2 ) = ^fc= 2 ( E ( n A:(9) 2 ) - E(n fc _i(< 2 ) 2 )) + E(IIi(g) 2 ) 
would diverge (logarithmically) as n —> oo. Thus if a limit exists it must 

7T 

equal M 2 ((/) = T(p q ) 1 ~‘ 2 T( q ) Jo 4 sin 29-2 0 cos 2q ~ 2 Odd (recall the notation of 


(59)). 


Let us thus write 


(61) E(n n (g) 2 ) = M 2 (g)-4 2 ). 

A priori, we only know that eh is a bounded sequence (as we showed in 
the previous section that E(II n (g) 2 ) is bounded for ^ < q < g 2 ), but let us 
( 2 ) 

try to show that eh —> 0. Plugging this ansatz into the recursion, one finds 

4+1 = 4 2) (l - h -(T{2q) - 2T(q))'j + S n , 
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where 5 n = 0((-) a ). We can of course write the solution to this recursion 

’’explicitly”: if we define So = e) ; and interpret the empty product (k = 
n — 1) as 1, then 


(62) 


n— 1 


n—1 


M - 


k =0 j=k -\-1 


II (l--m2q)-2T(q)) 


The product we can estimate by 


(63) 


n—1 


1- -(T{2q)-2T{q)) 

j=k +1 ' •’ 


i\ b(T(2q)—2T(q)) 

<c(-\ 

n ) 


where C is independent of k and n. Note that if b is large (but independent 
of n ), the terms in the product may be negative for small j, but this still 
is a valid asymptotic estimate for large n. Recalling that |e)fc| < Ck~ a , for 
some a > 1, we see that 


(64) 


n —2 

e (2)| < Cn -a + Cn -b(T{2q)-2T{q)) ^ k -<*+b(T(2q)-2T(q)) 

k =1 


< < 


Cn -b(T(2q)-2T{q)) ^ 
C n - b ( T ( 2 <i)- 2T ( ( i)) log 


n, 


Cn 1 ~ 


-a + b{T{2q)-2T{q)) < -1 
-a + b(T(2q)-2T(q)) = -l . 
-a + b{T(2q)-2T{q)) > -1 


( 2 ) 

Thus we see that for any fixed value of b, lim,,-^ Cn = 0 implying that 
E(II n (g) 2 ) M 2 (q). 


Let us now try to do something similar for E(II n (( 7 ) fc ). From the recursion 
for U n (q), one finds 


(E(U n+1 (q) k ) -E(U n (q) k )^ 


-( T(qk ) - kT(q))E(U n (q) k ) + ]T ^E(n n (^)E(n n (g) fc 

3 =1 K ' J ' 

i n do 


-3} 


x - / 9 -— sin 2 « 6 cos 2q{ - k ~ j) 6 + 0 

2 J 0 sin 2 6 cos 2 6 


a— 1 N 


With a similar argument as in the k = 2 case, we see that if E(n n (g) J ) 
were to converge for j < k, then the only possibility is to Mj(q ) which 
satisfies the recursion 


(65) 


M k (q) = 


1 


T(qk) - kT(q ) 


7 = 1 


i n 

X - 


3= 

d6 


2 J 0 sin 2 6 cos 2 6 


s\n 2qj 6 cos 2q[k ~ j) 6 
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with M\(q) = 1. 

Based on our experience in the k = 2 case, let us write 


( 66 ) 


K(U n (qy) = Mj(q) - e. 


(J) 


i,(j) | 


< CjUl a i 


and proceed inductively: let us assume that for j < k, 
for some aj > 0 independent of m. Plugging this into the recursion for 
E(II n (g) fc ) (and making use of Lemma 18 and Lemma 21), we find 


(67) e% = eW (l - ^(T(qk) - kT(g))^ + 0(n~ a ) 

for some a > 1 . Arguing as for k = 2, we see that this implies that \6n \ < 
Cn~ lk for some 7 *, > 0. We conclude that E(II n (g) fc ) -A- M k (q). 

□ 

6. Convergence of n n (g) for | < q < 1 

One can check that for ^ < q < 1, T(jpq) — pT(q) > 0 for any p > 0 
so by the previous section, E(II n (g) fc ) converges to M^{q) for all positive 
integers k. The natural question to ask then is does this imply that II n (g) 
converges in law. This is the case if the moments M^(g) uniquely determine a 
probability measure. This in turn can be checked by Carleman’s condition, 
or the following slightly weaker condition (for details, see e.g. Theorem 
3.3.12 and the discussion around it in m 

Proposition 29. Suppose the moments of the random variables X n satisfy 
E(-X^) —> fik cls n —> 00 for each k and 

1 J_ 

( 68 ) liiri sup — < 00 . 

k —^00 

Then X n converges in law to the unique distribution with moments pk- 

To check this property for the quantities Mk(q), we prove the following 
proposition: 

Proposition 30. There exists some C(q) >0 (independent of k) such that 
for\<q< 1, 


(69) M k (q) < ( C\q)) k k \. 

Proof. Let us make the ansatz Mj(q) < C J j\ for j < k. We then have by 
the recursion for M^ 


M k (q) < 


C k k\ 1 
T{qk) - kT(q) 2 

k 


< C k k\ 


T{qk) - kT(q) 2 


[ 4 cos 2 ^- 4 e V tan 2 ^'- 2 Ode 

1: 


cos 2qk ~ 4 9tan 2q ~ 2 9d0. 


Recalling that 
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(70) 


T(q) 


v^r(g-i) 

2 r(g-l)’ 


we see (e.g. by Stirling’s approximation) that T(qk ) ~ \fk for large k so 
that k/(T(qk) — kT(q )) ~ — 1/T(g) > 0 for q £ (1/2,1). Moreover, by say 
the monotone or dominated convergence theorem, the integral tends to zero 
as k —> oo. Thus there exists a fco such that for k > ko, 


<71) T( q k) - Wm So < 1. 

Also let us define 

(72) C7(g) = max f 3 

3<k 0 \ Jl J 

so that Mj(q ) < C(q)ij\ for j < k^. Our argument and choice of ko then 
implies that M k (q) < C(q) k k\ for all k. □ 

Making use of Stirling’s approximation, it is easy to check that this esti¬ 
mate is enough for us to be able to apply Proposition [29| and we get: 

Theorem 31. For 1 < q < 1, and any b > 0, n„(g) converges in law to 
the unique non-negative random variable whose positive integer moments are 
given by (M k (q)) k . 

Remark 32. As the moments M k are independent of b and the random 
variable whose moments these are is uniquely determined by them, the ran¬ 
dom variable is independent of b, so asymptotically, the only b dependence 
in P n (q ) is deterministic. 

7. The stationary equation for q e ( 5 , 1 ) 

In addition to describing the positive integer moments, one can also hope 
to describe for example the Laplace transform of the limiting law. Let us 
thus define <p n (t) = E(e“* n ”^) and let 4>(t) be the Laplace transform of the 
law of the limit (i.e. <f(t) = lim n _ ? . oc (j) n {t )). The main goal of this section is 
to prove the following proposition 

Proposition 33. For ^ < q < 1, cp(t) is the unique solution to the equation 

(73) T{q)t4>'{t) + ]- [ . 2 -— (c/)(t sm 2q d)4>(tcos 2q 6) - cj>(t)) = 0 

2 J o sin 6 cos z 0 

in the space of Laplace transforms of laws of non-negative random variables 
with unit mean. 

Remark 34. Note that our bound on the growth rate of the moments im¬ 
plies that (j) is analytic in some domain around the origin and this integral 
equation can be obtained from the recursion for the moments in this region. 
Nonetheless, we shall provide a proof valid for all t > 0. 
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Proof. Uniqueness follows from the fact that this integral equation implies 
that the positive integer moments satisfy the recursion of the previous 
sections and we saw that these moments determine the random variable 
uniquely. 

Let us write the recursion in terms of (j) n in the following way: 


(74) (fn+lit) - (fn{t) = E ( V™ (fn (t ™ ~ 0„(t)^ . 

Let us expand the integrand here: 


/ sin 2q 9\ ( cos 2q 9\ 


(75) 


\ — 10 \ .I / / 

= fn(t sin 2q 9)4> n (t cos 2q 9) - 4> n (t) 

+0 ”(‘lSy)^ (t ”- c)i ” s289 (sh) 

4- (h„(tc,os 2q 9)6'(t, 


- 1 




for some f n>c between the points tcos 2q 6 and and t, n _ s between the 

points tsm 2q 9 and 

We then note that the proof of Lemma [2] can be applied to functions 
depending on n as well as long as we have uniform bounds on the function 
as well as the uniform estimate f(9) = 0(\9\ a ) as 9 —> 0. Now it is easy to 
check that on [—7 t/4, 7t/ 4] one has 


(76) \4> n (t cos 2q 9)<f> n (t sm 2q 9) — (j> n (t)\ < tsm 2q 9 + t( 1 — cos 2? 9) 
so an application of Lemma [2] gives 

E(4> n (tcos 2q 9 n )(f n (t sin 2q 9 n ) - 4> n (t)) 

(77) = 2— f . (4>n(t cos 2q 9)4> n (t sin 29 9) - <t> n {t)) + o( (-*) 

n Jo sm 29 VW ) 

for a suitable a > 1. Note that by the uniform bounds on the integrand, 
the integral converges to the corresponding one with <f n replaced by 4>. 

For the other two terms, consider the function 

f(9) = -tT{q)4>' n {t) + 4> n (t (!>' n {t n)C )tcos 2q 9 

+ (j) n (t cos 2q 9)cf' n (t njS )tsm 2q 9. 

Noting that as (f>" a is uniformly bounded (as q < 1), we see that as 9 —> 0, 
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(78) f(6) = 0{\0\™ n{2 ' 2q) ) + O 

where both bounds are uniform in n. Thus making use of Lemma [2] once 
again, 

( t & r(j ) ^n(tn,c)t COS 2g 9 + (j) n (t COS 2g d)(l) , n (tn, s )tsm 2q B 

(79) = tT(g)<(t) + oQ), 
for some a > 1. We conclude that 



71 

-(0 n+ l(f) - <f> n (t)) 

— 1 f\ 

( 80 ) =\ [ . 9 A 4>n(tcos 2q 6) (j) n {tsm 2q 9) - <l> n (t))+tT{q)(l>' n (t) 

2 J o sin 26 



for some a > 1. Moreover, the bound is uniform in n. 


Now the left side converges as n —> oo (the integral by the dominated 
convergence theorem and the bounds we pointed out, and the derivative 
by the dominated convergence theorem when we write it as — IE(II n e~ fn ")). 
The left side can only converge to zero as otherwise 4> n would diverge log¬ 
arithmically, which is impossible. Thus we see that 0 satisfies the equation 
in the claim. 

□ 


8 . Discussion and open problems 

In this section we point out some issues we were not able to address and 
discuss some possible further questions related to the problem studied in 
this paper. 

8.1. Connection to band matrices. In Section [2] (see Remark|8| we dis¬ 
cussed a possible approach to going beyond the LME approximation for the 
band matrix model. A natural question would be to find an actual ran¬ 
dom matrix model for which a LME-type recursion becomes exact. For the 
analogous problem of Gaussian multiplicative chaos measures, i.e. measures 
whose density is given by exponential of the Gaussian Free Field the role 
of the LEM recursion is played by the Mandelbrot cascades resulting from 
replacing the GFF by a hierarchical version thereof. Note that the critical 
random matrix models with hierarchical structure do exist (26] 00], but that 
structure itself is not enough to make the LME recursion exact. Is there a 
hierarchical random matrix that produces a LME recursion? 
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8.2. Convergence for 1 < q < q c . Based on the analogy with the branch¬ 
ing random walk and the bounds on the moments E(n n (g) p ) for some p > 1, 
it is natural to expect that converges also in this regime. For q < 1, 

our approach centered around the fact that all moments existed and we were 
able to describe their limiting values through a recursion. 

For q > 1, all moments will no longer exist, but the stationary equation 
(73) still makes sense. For branching Brownian motion Bramson [12] (see 
also mi ei for the case of branching random walks) proved convergence by 
perturbing around stationary solution to the stationary equation. However, 
in that case nronotonicity due to maximum principle was an important in¬ 
gredient and we have not been able to prove such monotonicity in the LME 
case. 


8.3. Properties of the limiting random variable. In the case of the 
branching random walk a lot is known about the limiting random variable: 
it is almost surely positive, it has (in the Gaussian case) all negative moments 
and the tail of its distribution is known. In our case, our moment estimates 
suggest that the tail of the distribution might be similar i.e. we conjecture 
that for \ < q < q c the limit 

(81) lim x p ^^P(n(g) > x ) 

x—>oc 

exists. Moreover, the stationary equation suggests that if one writes n(g) 
for the limiting random variable, one should have P(n(( 7 ) = 0) = </>(oo) = 0 
and in analogy with BRW it is natural to conjecture n(g) almost surely 
positive and has all negative moments. 

8.4. q > q c and the freezing transition. The analogy with branching 
random walk leads one to ask whether T\ n (q c ) converges to zero and if so 
whether a nontrivial limit can be obtained with additional normalization. In 
the study of the branching random walk, Ii n {q c ) converges to zero under the 
case of a ’’second order phase transition” or the so-called ’’boundary case” 
(see e.g. M) which in our model should correspond to q c T'(q c ) — T(q c ) = 
0. There are branching random walks where this is not satisfied (though the 
critical point is characterized slightly differently) and their behavior for q > 
q c is significantly different from the ’’Gaussian case” see e.g. 0S|. In this 
case, n n {qc) converges to a non-trivial random variable without any further 
normalization. Thus for q c , the central question is does n n (g c ) converge to 
a non-trivial random variable or is the correct quantity to study — d q Ti n (q c ) 
(which can be defined e.g. using the tree construction). 

For the branching random walk the case corresponding to q > q c is con¬ 
trolled by the behavior of the maximum of the walk. In our case, using the 
tree construction, we can define a branching random walk which however 
is inhomogeneous and also the branches repulse each other. This leads to 
peculiar behavior. Indeed, consider the law of (log sin 2 9 n , log cos 2 9 n ). With 
probability of order 1 — ^, 9 n is of order 4 so log sin 2 6 is of order — log n 
while log cos 2 9 is of order —With probability of order ] n . 9 n is of order 
one so that both log sin 2 9 n and log cos 2 9 n are also of order one. This picture 
where apart from a rare event, the path of the maximum is deterministic is 
rather different from the standard branching random walk picture. While 
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studying this maximum is a non-trivial question, we refer to [24] , where it 
is argued that the scaling of the maximum is essentially the same as for a 
branching random walk with independent Gaussian increments. 

Finally, it is natural to ask is there a freezing transition phenomenon in 
the model. If there is a normalization under which II n (g) converges to say 
II(q) for q > q c , is this given in terms of the relevant critical random variable 
and a stable Levy subordinator? 

8.5. Geometry of the eigenvectors. The eigenvectors of the random ma¬ 
trix are given by 

'l>i r = O n i$, i = 0,...N-l 

where ('ipf)j = Sij • As in the derivation of the recursion © the recursion 
© on the resonant set can be used to arrive to the following recursion for 
the eigenvector. Let be the angles determined by the matrices AfT of 
(|8]) with j = i ± (n + 1). Fixing i either {i, i + n + 1} or {i,i — n — 1} 
belongs to the resonant set R n , with equal probability. Let cr” £ {0,1} be a 
Bernouilli (J, §-) random variable. The model for the eigenvector recursion 
is 

€ +1 = ( a i COS 9 n,i + ( 1 ~ a i) COS 

+<r? sin On^i+n+l + (! - a ?) Sin0"^* n -n-i 

This leads to the representation 

N -1 

(82) l/)f = ^ II Pn(u n ,u n+ 1) 

UJ 71=0 

where the sum is over all paths oj = (wo,wi, • • •, w^v) on Zv with ujq = i, 
UN = j and oj n+ 1 £ {w n , ui n + n + 1, uj n — n — l}and 

( aZ n cos e+ Un + (1 - O-SJ COS e~ Un Uj n+ 1 = u n 

Pn(u n , U} n+ 1 ) = l cr™ n sin 0+ Wn u n +i =u n + n+ 1 

( sin e n,u n w n+1 = u n - n - 1 

and 0 n .i. i £ Zv are i.i.d. copies of 6 n . 

It remains a challenge to use such a description for a more detailed under¬ 
standing of the structure of eigenvectors. 
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